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Abstract 

We connect and generalize Matiyasevich's identity #0102 with Bernoulli 
numbers and an identity of Candelpergher, Coppo and Delabaere on 
Ramanujan summation of the divergent series of the infinite sum of the 
harmonic numbers. The formulae are analytic continuation of Euler 
sums and lead to new recursion relations for derivatives of Bernoulli 
numbers. The techniques used are contour integration, generating 
functions and divergent series. 



1 Introduction 

In his personal journal on the Internet Matiyasevich ^1 has presented the 
identity #0102 with Bernoulli numbers: for even n > 2, 

E E @ (i) 

k+l=n k+l=n \ / 

and welcomed readers of his website to volunteer further information on the 
topic. Here, B n J s are the Bernoulli numbers given by 



= Y,(- l T B n- (2) 
1 ^— 1 n\ 



n=0 

and H n = Sj=i j is the n th harmonic number. Note that the definition of 
Bernoulli numbers given above [2H] changes only the value of B\ from the 



usual definition and Matiyasevich's identity is not affected by this modifica- 
tion. In this paper, we use the well known tools of generating functions 22 , 
contour integration and divergent series to connect Matiyasevich's dis- 
covery to two other classes of identities. The first one due to Candelpergher, 
Coppo and Delabaere jH] is given by 

m 3 11 

£ H n = ty + - log(2vr) , (3) 

n>l 

where [7Z] represents the Ramanujan summation and 7 is Euler-Mascheroni 
constant. 

The other class of identities related to Matiyasevich's discovery is given 
by explicit formulae for Euler sums in terms of Riemann-zeta values. We 
give a sample of typical examples of such identities below. 

OO TT 

= 2C0) (4) 

n=l 

£§ - \m w 

n=l 

OO jT 

= 3C(5) - C(2)C(3) (6) 

n=l 

The left hand side of these identities are examples of Euler sums. For the 
definition of general class of Euler sums and their formulae, see ^U]. Apostol 
and Vu [21 and Matsuoka ^S] studied the related Dirichlet series h(s) = 
T,n=i ffS > 1 and showed that h(s) can be analytically continued to 
the entire complex plane. For the history of Euler sums, see [3] and for a 
recent survey see jl] and the references therein. Proof and generalizations 
of Matiyasevich's identity and its relation to Miki's identity are discussed in 
H3, d, QH and 0. 



2 A Hybrid of Matiyasevich and Candelpergher- 
Coppo-Delabaere Identities 

In this section we prove the following result which combines the results of 
Matiyasevich and Candelpergher-Coppo-Delabaere. 
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For n > 1, 



h(-n+l) = B' n + nB n _ 1 + 7 B n - ]T ) ^B, - (7) 

k+l=n \ I 



where h(—k) = Ylr^=i H n n k is understood in the sense of an appropriate 
sum of a divergent series. 

Remarks 1. The last two terms are identical to those in Matiyasevich's 
identity (JTJ. 

2. The first identity in this set is that of Canderpergher,Coppo and De- 

1 l 

labaere. By substituting n = 1 and using the fact [23] that B[ = — — — log 2ir 

we get 

oo 

h(0) = J2H n = -~f + ---log(2n). (8) 

n=l 

The difference between the identities (J3J and (JHJ) comes from the fact that 

m 1 

the Ramanujan sum of the divergent series 1 = — whereas the value 

n>l 2 

arrived at by standard summation methods is — . Such differences often 

J 2 
happen when differing methods of summing divergent series are used. See 

for example Chapter 5 of 20 for a recent exposition. 

Proof. Let us take the term 

in Matiyasevich's identity (JTJ). Recall that [22], if a(z) = ^ ( ^^- z n and 

^ b(n) . \ c ( n ) n i c ( n ) a (k) b{l) 

b ( z )=Y,^T z n ,thena(z)b(z) = j2^rz ™, where = £ -^T^p 

fc+Z=?i 

Choosing the generating functions 

z n 



1 ^— 1 n! 



n=0 

(set z — » — z in ([2])) and 

— L ) = E — £ T. (ii) 

z / , n n! 

n=l 
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see [T2] and [3], one gets ©. That is, 

z z - 1 



v v " '•• 

n=l fc+i=r 



k k\ l\ 



log 



(log (e 2 - 1) - log z) 

(-,) 1 ° g(e '- 1) + - 
e z — 1 e * 

( ^ loge z (l-e~ z ) + 



log 2 



log(l — e z ) 
1 - e~ z 



+ 



1 

ze 



log z 

z 

— log 2- 



(12) 



In order to pick out the n th coefficient of the series in (|12j) . let us integrate 
along the Hankel contour C. C is a loop that starts at infinity on the 
negative real axis, encircles the origin in the positive direction, excluding 
the points ±2z7r, ±4i7r • • •, and returns to negative infinity. Thus, 



E 



Bh B, 



1 



k k\ 11 



k-\-l=n 

and from (|T2*|). we get 

Bk Bj 
k k\ 11 



1 



Bk B, 



2m Jc z n+1 ^ , ^ k k\ /! 

n=l k+l=n 



z n dz . 



(13) 



E 

k+l=n 



1 



1 



27ri Jc z n e z — 1 

1 log(l - e~ 2 ) 



+ 



+ 



1 



27TZ JC Z n 1 



1 



1 



logz dz 



2m Jc z n 1 — e 

£ n _! (_i)«-ifc(_ n + i) a 



(n- 1)! 



(n-1)! 



(14) 



(15) 



where -B^ = ^_B s | s=n is the derivative of the Bernoulli number. Multiplying 
both sides of ()15|) by n! gives (J7J) which is an extension of Matiyasevich's 
identity #0102. The evaluation of the three integrals is given below as 
Lemmas 1, 2 and 3. 
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Lemma 1 For n > 1, 



-/-^^ = AV (16) 
2m Jc z n e~ z - 1 {n - 1)! 

Proof. Immediate from (flOJl . 

Lemma 2 For n > 1 

1 f 1 log(l - e" 2 ) ^ = (-l) n h(-n+l) 



2m Jc z n 1 — e z (n — 1)! 

Proof The integrals along the upper and lower lips of the Hankel con- 
tour C cancel each other and thus the integral is equal to the residue of 
1 log(l - e~ z ) 

at the origin. We now evaluate this residue. Using the 

z n 1 — e~ z 
identity 

-!5f^=£ff„x», (18) 

1 — X , 

n=l 

where H n is the n th harmonic number defined earlier, we get 

log(l — e" 



E H n e- nz . (19) 



1 - e- z 

n=l 

Expanding the exponentials as Taylor series around the origin, and rear- 
ranging the divergent series we get 

E z n e-~ = E E Hn { -^- = E ^-f-K-k) • ( 2 °) 

n=l n=lfc=0 k=0 

Lemma 3 For n > 1 

1 r 1 p 2 R' R 

TW 4t7^ logz^ = ^-^(-7 + ^n). (21) 
Z7rz Jc z n 1 — e z n! n! 

Proof See Lemma 1 of [H]. We have P 

5 s+ i 1 /■ 1 



cfe . (22) 



T(s + 2) 2vri Jc z s+1 1 - e 
By taking derivatives on both sides of (|22|) with respect to s, we get 
B' s+1 B s+1 T'(s + 2) 1 y 1 e 



T(s + 2) T(s + 2) T(s + 2) 2m J c z s+1 1 - e z 
Using the identity [7] 



log z dz . (23) 



r(n + i) 

and putting s = n — 1 in (|23j) , we get (|21l) . 



r'(n+l) , s 

+ (24) 
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3 Possible Directions for Generalization 



Let us take the integral 



roo T s— 1 poo , . 

/ — dx = / X s " 1 ( e~ 2x + 2e~ 3x + 3e~ Ax + • • • ) x 3 ' 1 dx 

Jo [e x - I) 2 Jo v ' 

12 3 \ , . 

2^ + ^ + 4^ + + ---J r(s) 
2-1 3-1 4-1 \ . , 

— + — + — + + '-j r ^ 
= T(s) (C(s - 1) - COO) , (25) 

if 5?(s) > 2 |21j . Writing this real integral as complex integral using the 
Hankel contour gives 

C (,- 1 )- {W=r(1 - s) ^L| c |^ d2 (26) 

The integral on the right hand side can be analytically continued to the 
entire complex plane as is done in JJ|, and now replacing s by — s in 1)26(1 we 

C(-s-l)-g-s) _ J_ /■ 1 e 2 * 

r(l + a) 2iri Jc z s + l (1 - e z ) 2 [ ' 

Differentiating both sides with respect to s and writing log z = log(e z — 1) — 
log( e ~ 1 ), the right hand side of the integral at s = n will give, 

1 / 1 e 22 logz & 1 f 1 fe z logz\ / e* \ 



2vri z n+1 (1 - e z ) 2 ' 2vri ^ n+1 V e z - 1 / Ve z - 1 

27ri 7c ^ n+1 V e 21 - 1/ v e* - 1 / 

1 ' 1 w^v^v 

(28) 



2ni Jc z n+1 V V z J J \e z — 1 



Writing the power series for the two terms in terms of h(—k) and B n will 
give relationship between B' n+ ±, B' n and h(—k) for k < n. As we have 
already derived a connection between h(—n + 1) and B' n in Sectional this 
would lead to a recursion relation for B' n and Glaisher-Kinkelin-Bendersky 
constants |13j , jSj . It has not escaped our attention that finding a recursion 
for B' n will lead to a formula for C(2& + 1) from the functional equation of 
the ("-function. 
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4 Comments 



There is a lot of confusion in the literature regarding Riemann zeta 
functions, Bernoulli numbers and their representations. There are es- 
sentially real and complex representations of Bernoulli numbers and 
Riemann zeta functions. We are inspired by Rota to use the idea 
that Bernoulli numbers can be generated by a complex representation 
around the Hankel contour. The idea of using generating functions for 
Euler sums was powerfully demonstrated by Flajolet. Gessel's deriva- 
tion of Miki's identities is on similar lines but cancels the logarithmic 
part. Roman and Rota in JH] develop the analogue of logarithmic 
Taylor series and attempt to sum derivatives of Bernoulli numbers. 
This paper essentially combines the ideas of these authors. 

The Nor Kind polynomials Bn^ are defined by 



71=0 

Note that that Matiyasevich's identity (^Q) is related to the derivative 
of the above at z = 1. 

3. As Euler sums are related to zeta functions in several identities as 
an after thought it seems obvious that using the functional equations 
for Euler sums and zeta functions, there have to be several identities 
which are analytically continued versions of existing identities. 

4. In 9 Dunne and Schubert use the generating function 

1 00 B 1 

fc) = *W ~ lo % x + Tx ~ £ it j» ' (30) 

k=l 

where ip(x) = ■ If the logarithm term is not killed, using the 
ideas in ^0], new identities may be obtained by a new method. These 
identities should also be related to analytic continuation of Euler sums. 

5. In jH] the importance of the combination £'(— n) — H n (^(—n) in the the- 
ory of Glaisher-Kinkelin-Berndersky constants is stressed. This com- 
bination occurs naturally in this derivation. 
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6. A novel feature of the approach outlined is the surprising appearance 
of a divergent series as a residue. These formulae cannot be checked 
easily numerically as they involve divergent series. Further there seems 
to be some difference between the conclusions of Apostol and Vu [2] 
and Candelpergher, Coppo and Delabaere regarding the values of 
the analytically continued Euler sums at negative integers. This would 
require very rigorous and careful analysis. Hence the derivations we 
have given should be considered heuristic and merely sketching the 
rich possibilities existing in this direction. 
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